Skyrme crystal or Skyrme liquid? 
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We address the quantum melting phase transition of the Skyrme crystaL Based on generic sum 
rules for two-dimensional, isotropic electron quantum liquids in the lowest Landau level, we propose 
analytic expressions for the pair distribution functions of spin-polarized and spin-unpolarized liquid 
phases at filling factors 2/3 < i' < 1. From the pair distribution functions we calculate the energy 
of such liquid phases and compare with the energy of the solid phase. The comparison suggests 
that the quantum melting phase transition may lie much closer to !/ = 1 than ever expected. 

PACS: 73.40.Hm 
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Whether or not a solid phase composed of skyrmions 
[Q exists near filling factor ly = 1 poses, probably, one 
of the most fascinating problems in recent history of the 
quantum Hall effect [g[. This phase, known as Skyrme 
crystal |Q], is, nowadays, the target of a number of ex- 
periments based on various techniques which, in the past, 
have borne out the existence of the skyrmions themselves 
[Q-g[. To date, however, Uttle experimental evidence of 
the existence of the Skyrme crystal has been provided. 
Bayot et al. Q have reported a peak in the specific heat 
at very low temperatures which may be attributed to 
a phase transition from a finite-temperature liquid phase 
into a solid phase at zero temperature p,|7|| . The spin po- 
larization measured by Barret et al. Beseems to be con- 
sistent, at least in the range of filling factors |l"i^| ~ 0.2, 
with that calculated by Brey et al. |^ for a proposed 
solid phase. This apparent consistency neither guaran- 
tees that a zero-temperature liquid phase exists out of 
that range nor, most importantly, rules out its existence 
close to i^ = 1; for an alternative proposal of an unpolar- 
ized or partially polarized quantum liquid is not available 
to compare with. 

If a solid phase really exists at zero temperature, the 
following question needs to be answered: What is the fill- 
ing factor at which quantum fluctuations turn the solid 
into a liquid?. Filling factors i/ = 2/3 and v — 4/3 set 
the lower and upper melting transition limits, respec- 
tively. (We will focus on i^ < 1 since the physics for 
I' > 1 is similar due to particle- hole symmetry.) For 
large enough Zeeman energy per electron, gfiBB/2, the 
v = 2/3 state is a spin-polarized liquid. More precisely, 
it is a i^ = 1/3 Laughlin state of holes. For g^BB/2 — 
there has been mounting numerical and experimental ev- 
idence of a spin-unpolarized liquid ground state g]. To 
zero in on the actual location of the melting transition 
one can try to estimate the energy scale of the quantum 
fluctuations associated to the phonon modes in the clas- 
sical solid phase H] and compare it with some relevant 
energy scale, typically, set by the Coulomb interaction. 



This is a very crude approximation; moreover, it does 
not provide any information about the liquid phase. An- 
other approach, used in Wigner crystal studies |Q] and 
which we explore in this Letter, is to offer an alternative 
description of the ground state in terms of a quantum 
liquid, obtain its energy, and compare it with that of the 
solid phase. 

After Laughlin's proposal of a wave function for the 
main filling factors v = \/M (where M is an odd integer 
> 3), many other wave functions have been proposed for 
other filling factors with notorious success [Q]. For many 
purposes, however, a many-body wave function contains 
redundant information which can be compacted in a more 
advantageous way. In this spirit we note that the ground 
state energy per particle, E /N = e, for an infinite, neu- 
tral, and isotropic liquid state can be obtained from the 
pair distribution function g(r) [^: 



dr r[g(r) — l]y(r). 



(1) 



where V{r) is the interaction potential. In general, g{r) 
can only be obtained from the wave function using the 
plasma analogy pj. On occasions |lO|, however, one can 
find an analytic approximation for g{r) based on the fact 
that it must fulfill generic sum rules. In addition to 
these sum rules and in order to get the desired accuracy, 
one must provide microscopic information specific to the 
quantum liquid state. This information is essential, but, 
as shown below, can be obtained from simple consider- 
ations which do not necessarily involve the wave func- 
tion. Here we find g{r) and calculate the energy of spin- 
polarized and spin-unpolarized quantum liquid phases in 
the range of filling factors between v = 2/3 and v = \. 
Exact diagonalizations for v — 2/3, where we certainly 
expect liquid phases, and v = 4/5 allow us to confirm the 
reliability of this procedure. Figure y shows the energy 
per particle of both liquid phases compared to that of the 
solid phase proposed by Brey et al. in a mean-field ap- 
proximation. Our results confirm the existence of liquid 



phase ground states at i^ = 2/3 and also at v — 4/5 [ plj . 
Most importantly, our results suggest that, even without 
considering partial polarization of the liquid phase, the 
quantum melting phase transition takes place very close 
to :^ = 1. In fact, it seems to occur near the transition 
point that separates the square Skyrme lattice from the 
triangular one 
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FIG. 1. Energy per particle in the Hartree-Fock approx- 
imation of the solid phase or Skyrme crystal (dots) and 
of the polarized (stars) and unpolarized (solid line) liquid 
phases proposed in the text. For completeness, the values 
obtained from the exact diagonalizations are presented by cir- 
cles for both polarized and unpolarized phases. The results 
for the polarized liquid phase include a Zeeman energy shift 
oi qHbB /2 — — 0.0075e^/£ and all the energies are referred to 

Following early work by Girvin [|0| , we start by noting 
that the pair distribution function of a liquid state in the 
lowest Landau level (LLL) can be expanded in terms of 
wave functions for relative states of two particles. Girvin 
noticed that it is convenient to recast this expansion in 
terms of the pair distribution function for the i/ = 1 
liquid plus a perturbative series. When the liquid is fully 
polarized 



m=l ' ^ 



(2) 



where the prime denotes a sum restricted to odd val- 
ues of m. The generalized g{r) for fully spin-unpolarized 
states (i.e., with equal density of spin- up and spin-down 
electrons) 



gif) = t^[.9TtW +.9TiW]. 



(3) 



is divided into a contribution coming from equal spin 
electrons, g||(r)[= .9xx (»")], and a contribution from op- 
posite spin electrons, 5fx(r)[= ^^(r)]. It is natural to 
expand g^^ir) in a way similar to that in (pi). It is also 



convenient to make a similar expansion for 5|x(r), ex- 
cept for the fact that the sum runs now over all possible 
non-negative integer values of m: 



9-\i{r) 



1 



-r^/2 



-r^jA 



m=0 



(4) 



If, for all TO, Cm ^ and c„i —* 0, one recovers the 
partial distribution functions 5x|(r)/2 and g-\i[r)l2 of 
the V = \ liquid ferromagnet with minimum polarization 
in the direction of the field (S'z = 0). Note that, in 
general, not any set of coefhcients Cm and c,„ represents 
a good liquid. First of all, at any filling factor, c™ — > 
and Cm — > for to — > oo since g(r) —^ 1 for r — > 
oo. Secondly, the pair distribution function must satisfy 
certain sum rules which establish relations between the 
coefficients. Finally, these coefficients must be consistent 
with the microscopic information one must provide. 

Spin-polarized liquids. When one applies the charge 
neutrality, perfect screening, and compressibility sum 
rules IQ to the expression (g) one obtains the follow- 
ing set of equations: 



J2 'c™-(l-^-^)/4 

m— 1 

oc 

J2 '(m-t-l)cm = (l-i/-i)/8 

m— 1 

oo 

Y^ '(TO-t-l)(TO-|-2)Cm = (l-i.-l) 



(5) 
(6) 
(7) 



The charge neutrality and perfect screening sum rules 
can actually be derived exactly from general considera- 
tions applicable to any isotropic liquid state in the LLL 
|l2[ . In short, they come about due to the fact that the 
number of particles and the total angular momentum are 
good quantum numbers for the Hamiltonian. The com- 
pressibility sum rule can only be invoked as it is if there is 
a classical plasma analog similar to that for the Laugh- 
ling states M. For the filling factors considered below 
this analog exists Q . 

Now, in order to incorporate the microscopic informa- 
tion related to the short-distance behavior into the pair 
distribution function, we first express it in a more ortho- 
dox way [O : 



9{r) = — e"'' ^^ ^ {rimno) 



2™to! ' 



(8) 



A series expansion of (|^) allows us to relate the coeffi- 
cients Cm to the correlation factors {nmno) = (cJ„CmcJco) 
which we try to calculate in what follows. The particle- 
hole symmetry dual filling factors of z^ = ^/M, v = 
1 — l/M , lie in the range we want to explore here and we 
assume a Laughlin state of holes at these filling factors. 



There are no simple wave functions representing these 
Uquids, but the same particle- hole symmetry allows us 
to write 

{rimno) = 1 - {hl^hra) - {hlho) + {h\^h„-,hlho) (9) 

where the h\^ and hm represent creation and annihilation 
operators of holes. For the proposed Laughlin states 



(hl^hmhlho) =0, m<M-l 



(10) 



which gives {nmno) = 2;^ — 1 for to < M — 1. With this 
result we find Cm — —(1 — i^^^)^ for odd to, < A/ — 2 
and, with Eqs. rtq), (0), and (R), we calculate the next 
three coefficients, setting the rest to zero. In Table ffl 
we present the energy calculated from the analytic g{rys 
so obtained for z/ = 2/3 and ly = 4/5 along with that 
obtained from exact numerical diagonalizations on the 
spherical geometry for the same filling factors. The per- 
fect agreement between these two values confirms the ex- 
istence of a Laughlin liquid state of holes at those filling 
factors and, moreover, confirms the reliability of the an- 
alytic pair distribution functions. We calculate now the 
energy of Laughlin liquids at any filling factor beyond 
i> = 4/5 within the set ly = 1 — 1/M (stars in Fig. |^). 
As t^ ^ 1, numerics cannot be employed as a test of 
accuracy. However, all the coefficients must vanish in 
this limit and we expect the expansion (0) to get even 
more reliable. In fact, our results are consistent with the 
known behavior e{i' = 1) = de{i')/di'\j^^i- ~ — i-\/7r/2 
PI and with the expected fact that the triangular Skyrme 
crystal (the two closest points to i^ = 1 in the sofid phase) 
has lower energy than the Laughlin-like liquid states near 
V ~ 1. Figured also shows that the polarized liquid phase 
beats the solid phase aX v = 2/3, but not dX v = 4/5, in 
agreement with Ref. p]. 

Spin-unpolarized liquids. The two first previously men- 
tioned sum rules can also be invoked for the pair distribu- 
tion function (^) of a two-component liquid. We obtain 
three sets of equations: 



^ c„ = (1 - 2i.-i)/4, ^£„ = l/4, (11) 

171=0 m— 

^(m + l)(c„ + c„) = (l-i.-i)/4, (12) 

where we implicitly consider Cm = for to, = 0,2,4,6,... 
The first two equations refer to g^i and g^i^ separately. 
They result again from generic considerations p^ which 
take into account the fact that the number of particles 
for each component (spin orientation) is a good quantum 
number. These two can be combined into a single one 
since the energy of the unpolarized phase depends only 
on 5|| -I- (7||, i.e., on c^ -I- c^. The compressibility sum 
rule has not been invoked this time for we are not aware of 
any two-especies plasma analog for the filling factors con- 
sidered above. Its existence cannot be discarded though. 



Next, we provide the microscopic details of the unpo- 
larized liquid states we propose. An expansion like (|^) 
including the spin degree of freedom allows us again to 
relate the coefficients Cm and Cm to the correlation factors 
("mT'^OT)(= {nTninoi)) and (n™|noi)(= ('^minot))- The 
duality between v and 1 — z^ in previous section allowed 
us to propose sensible liquid ground state wave functions 
at filling factors i/ = 1 — 1/M, but, here, the spin de- 
gree of freedom breaks such a duality. However, note 
that one can always find an appropriate hard-core inter- 
action model for which the polarized ground state satis- 
fies Eq. {Kw, whether or not its wave function is explic- 
itly known. Similarly, the hard-core interaction model 
Vb^Vi>V2 = V3... = guarantees that the unpolar- 
ized ground state satisfies (n-ofnoi) — (n-oi^iot) = for 
V < 1. It also guarantees the upper limit 

(nino) = 2(ni|noT) + 2(ni|nox) < {2v - 1), (13) 

which, as follows from previous arguments, is the value 
associated to the fully polarized state at filling factors 
2/3 < V <\ ioY any Sz- The inequality refiects the free- 
dom for a pair of particles to occupy states with even 
relative angular momenta. We propose the following fac- 
torization: (niTio) = f{v){2v — 1) with f{v) < 1. In 
general, /(z^) = a + (1 — a)v since we do not expect 
terms proportional to v^ or higher to be relevant in a 
density-density correlation function and the condition 
f{v — I) — 1 must be fulfilled. Thus we set cq = 
and ci + ci = 2(1 - a)(l - v-'^) - 2a{l - v'^f and, 
using Eqs. (|l^) and (|l2|), we can calculate C2 and C3 -t-ca. 
We especifically set a = 25/8 — 2\/2 which builds a g{r) 
from which we reproduce the expected behavior of the 
energy for a very dilute state of skyrmions at vanishing 
Zeemancoupling, de(i^)/dz^|^^i- = — |;-\/7r/2. Moreover, 
this choice of a is also consistent with our best numerical 
estimate of (nino) for v = 2/3. As Table | shows, the 
energy obtained from the analytical expression of g{r) 
for V = 2/3 and v = 4/5 and the exact results for the 
real Coulomb interaction are remarkably close. The solid 
phase energy (see Fig. ||) lies clearly above the exact and 
estimated values, which proves the existence of unpolar- 
ized liquid phases at these filling factors. (Note that the 
energy of the polarized phases are always above the corre- 
sponding ones of the unpolarized phases, but a partially 
polarized liquid at v — 4/5 is possible and might have 
lower energy for the value of the Zeeman coupling con- 
sidered). One can interpret these unpolarized states as 
ly = 1/3 and ly = 1/5 hole liquid states, with each hole 
being "accompanied" by one and two spin flips (i.e., spin 
waves), respectively. Although we cannot present a rig- 
orous argument, it is tempting to interpret this hole-spin 
waves association as a skyrmion nm , and call these states 
Skyrme liquids. 

In order to address the quantum melting phase transi- 
tion we need the energy at any filling factor 2/3 < ly < 1: 
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Although, as required, this expression reproduces known 
results fairly well, we have to stress that its derivative is 
continuous as a function of v. In other words, it does not 
capture the expected gaps in the chemical potential at 
the main filling fractions. Experimentally, however, only 
the gap dX V = 2/3 is clearly visible for fully-polarized 
phases. Towards v — 1 the gaps &i v = 5/7 and v = 4/5 
are the next in importance, but they barely show up in 
the longitudinal or Hall resistivities even at the lowest 
temperatures. In addition, the gaps are smaller in unpo- 
larized phases B so we expect the analytical expression 
(O) to capture the overall behavior of the real e(i^) of 
a spin-unpolarized liquid, particularly close to z^ = 1. 
Before drawing conclusions about the quantum melting 
phase transition from Fig. n^one should keep in mind that 
quantum fluctuations have not been considered in the 
solid phase calculation and that they are bound to lower 
the energy of this phase. The magnitude of this correc- 
tion has been carefully estimated for the Wigner crystal 
and it is not unreasonable to expect an energy shift of 
the same order of magnitude for the Skyrme crystal. Bor- 
rowing the results from Ref. Q we estimate the shift to 
be — 0.001e^/£ for v = 4/5, decreasing in absolute value 
as V -^ \. In our favor we must also say that partially 
polarized liquid phases, which have not been considered 
here, are expected to lower the ground state energy for 
typical values of the Zeeman energy. All considered, our 
results seem to place the quantum melting phase tran- 
sition very close to ;/ = 1. As expected, the triangular 
Skyrme crystal dominates sufficiently close to :/ = 1, but 
the transition point to the square lattice {v « 0.96) lies 
near the one where the unpolarized liquid phase begins 
to take over [u w 0.92). 

There are still many open questions regarding the na- 
ture of the zero-temperature melting phase transition 
over which this work does not shed any light. Disorder, 
for instance, plays a major role in this matter and has 
not been considered. Still, our results cast a caveat in 
what refers to the range of existence of the square-lattice 
Skyrme crystal and, in addition, provide an alternative 
description for the electronic ground states in the range 
of filling factors 2/3 <v<l. 
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V 


iV = 4 Ar = 6 Ar = 8Af = 10Af = 12Ar-.cx) 


5(r) 


2/3 (P) 
2/3 (U) 


-0.5501 -0.5396 -0.5341 -0.5310 -0.5289 -0.518 
-0.6010 -0.5744 -0.5618 -0.5545 -0.527 


-0.518 
-0.528 


4/5 (P) 
4/5 (U) 


-0.5804 -0.5649 -0.5607 -0.552 
-0.6580 -0.6069 -0.5923 -0.56(7) 


-0.552 
-0.566 



TABLE I. Energy per particle (in units of e^/i") obtained 
from exact diagonalizations on a sphere for finite number of 
particles, A'^, at filling factors v = 2/3 and v — 4/5 (when 
possible, values for Af > 12 have been obtained, but are not 
shown here). The spin-polarized and spin-unpolarized calcu- 
lations are denoted by P and U, respectively. We also include 
the thermodynamic value obtained from a 1/A'^ cuadratic ex- 
trapolation and the energy obtained from the pair distribution 
functions described in the text. The last digit of the extrapo- 
lated value for v = 4/5 (U) cannot be taken at face value due 
to the small number of points available for the extrapolation. 



